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Abstract. Evolution of the cluster population has been 
recognized as a powerful cosmological tool. While the pre- 
sent-day abundance of X-ray clusters is degenerate in 
as, n and Qq, Oukbir and Blanchard (1992, 1997) have 
pointed out that the number density evolution of X-ray 
clusters with redshift can be used to determine Sl - Here, 
we clarify the origin of this statement by identifying those 
parameters to which the evolution of cluster number den- 
sity is most sensitive. We find that the evolution is con- 
trolled by only two parameters: the amplitude of fluctu- 
ations, crjvf) on the scale associated with the mass under 
consideration, R — 9.5/i 1 / 3 f2 ( ^ 3 M^ 3 /i _1 Mpc, and the 
cosmological background density, Q . In contrast, evolu- 
tion is remarkably insensitive to the slope of the power 
spectrum. We verify that the number density evolution of 
clusters is a powerful probe of the mean density of the uni- 
verse, under the condition that <tm is chosen to reproduce 
current-day abundances. Comparison of the cluster abun- 
dance at z ~ 0.5 — 0.6, from the EMSS, to the present-day 
abundance, from the ROSAT BCS sample, unambiguously 
reveals the existence of significant negative evolution. This 
number evolution, in conjunction with the absence of any 
negative evolution in the luminosity-temperature relation, 
provides robust evidence in favor of a critical density uni- 
verse (Oo = 1), in agreement with the analysis by Sadat 
et al. (1998). 

Key words: Cosmology: observations - Cosmology: the- 
ory - large-scale structure of the Universe - Galaxies: clus- 
ters: general 



1. Introduction 

X-ray galaxy clusters offer several interesting ways to con- 
strain cosmological parameters. The temperature of the 
intra-cluster gas can be related to the virial mass accord- 
ing to 



T = 4Af 1 2 5 /3 (l + z)keV 
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where Mi 5 = M/10 15 M© and we hereafter assume h — 
i?o/100.kms _1 Mpc _1 — 0.5. Such a relation can be easily 
deduced from the equation of hydrostatic equilibrium for 
the gas, leading to a temperature some 20% larger than 
the value given in Eq. 1, which was inferred from numeri- 
cal simulations (Evrard, Metzler & Navarro 1996). Typical 
clusters have temperatures between 2 and 14 keV, corre- 
sponding to scales between 5 and 15 h l / 5 VL ft, _1 Mpc. 
Henry and Arnaud (1991) have shown how both the nor- 
malization and the slope of the power spectrum can be 
inferred from the local temperature distribution function, 
a technique which has been widely employed in recent 
years (see Bartlett 1997 for a review). Oukbir & Blan- 
chard (1992) proposed that the evolution of the X-ray 
cluster abundance is a powerful probe of the mean cosmo- 
logical density, fl^. In order to apply the technique, Oukbir 
& Blanchard (1997, hereafter OB) established a detailed 
description of X-ray clusters. This new approach based on 
evolution has also received much attention lately (Henry 
1997; Bahcall et al. 1997); however, because the evolu- 
tion of cluster number density depends in principle on the 
mass considered, the spectrum of the initial fluctuations, 
its normalization and the cosmological framework, doubts 
have been raised concerning the validity of the technique 
(Colafrancesco et al. 1997). The purpose of this letter is to 
clearly identify the parameters controlling cluster number 
density evolution and to examine what one may say about 
fio using current data. 

2. Cluster properties and the mean density of the 
universe 

The Press-Schechter (1974) formalism, PS hereafter, is a 
rather simple description of the mass function and its evo- 
lution, and it has been shown to be in good agreement 
with numerical simulations (Lacey & Cole 1994). The PS 
prescription is: 



N{M, z)dM 



2 p c ilo S C A diner 
7T M 2 V^dhiM 



exp 



(S c Af 



dM 



where A = A(£lo,z) is the growth rate of linear den- 
sity perturbations, p c is the Einstein-de Sitter density, 
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Fig. 1. In this figure, we plot the evolution of the relative 
abundance of clusters, expressed by equation as a func- 
tion of redshift. The continuous lines are for M%5 — 0.125; 
the dashed lines, for M15 = 1.0; and the dot-dashed 
lines, for A/15 = 6.55. The thin lines correspond to an 
n = power spectrum index, while the thick lines are for 
n = — 2. The upper graph represents f^o = 1, and the bot- 
tom is for 17o = 0.2. In the first case, erg varies from 0.35cr m 
to 2.56er m ; in the second case, it changes from 0.82tr m to 
5.5<r m . Each figure illustrates that the relative abundance 
of clusters does not depend on the spectrum, nor on <7g. 
Clearly, the only important parameters are a m and the 
cosmological density parameter, f2 , the latter evinced by 
the difference between the top and bottom panels. 



S c = S c (flo,z) is the critical linear over-density required 
for collapse and <jm is the present-day amplitude of den- 
sity perturbations on a scale M. From this expression, it 
is clear that the cluster abundance at redshift z = 0, for 
a given mass M, determines <tm almost independently of 
the value of the spectral index, n, or of the density param- 
eter. There is only a slight dependence on these quantities 
due to their presence in the pre-factor of the exponen- 
tial term (and there is almost no influence from a possible 
cosmological constant). In practice, matching the present- 



day number of observed clusters with T > 4 keV requires 
a-(10 15 Af©) ~ 0.6 in an f^o = 1 universe, and a similar 
value, <t(10 15 M q ) ~ 0.8, in an f^o = 0.2 universe. How- 
ever, it should be kept in mind that this corresponds to 
two different linear scales of the initial density perturba- 
tion field: 



R w 9.5h 1/3 rtg 1/3 M^ h^Mpc 



(2) 



the difference being almost a factor of 2 between J7o = 
1 and Hq = 0.2. This means that the abundance of 
10 15 Mq clusters determines the amplitude on different lin- 
ear scales. In an Hq = 1 cosmology, a is fixed on a scale 
of 8/i~ 1 Mpc, the traditional normalization scale, while in 
an Sin = 0.2 cosmology, a is instead set on a scale of 
15/i _1 Mpc. Accordingly, for an open model (f2o ~ 0.2), 
cr(8/i Mpc) must be found by extrapolation using a spe- 
cific n, and is uncertain by a factor of two (see OB, Fig. 

1)- 

Let's now examine what governs the redshift evolution 
of the cluster abundance on a given mass scale, M. As in- 
spection of Eq. 1 clearly shows, only ctm, on the mass-scale 
considered, A and S c govern the evolution with redshift. 
As the the later two quantities only depend on f2, the 
redshift evolution is completely independent of the power 
spectrum index, n, and does not depend explicitly on the 
normalization at 8h~ 1 Mpc. This makes the redshift evo- 
lution remarkably simple to understand and to employ as 
a cosmological probe: once tr(M) is set by the present- 
day cluster abundance, the evolution of the number of 
clusters on the same mass scale is entirely and uniquely 
determined by the cosmological background (Oq). This is 
the essential reason for the robustness of the cosmological 
test originally proposed by Oukbir and Blanchard (1992). 
In order to illustrate this point, we define the quantity 



C m {z) = 



n(> M, z) 
i(> M,z = 0.) 



(3) 



as a simple measure of redshift evolution. We plot C m (z) 
in Fig. 1 for different spectra normalized to the same am- 
plitude, cr(M), and for two different cosmological back- 
ground densities - il = 1 and f2 = 0.2. It is important to 
note that, for fixed cr(-M), as varies as n and M change; 
for example, when o~m = 1, the normalization eg goes 
from 0.3 to 1.3 as n is varied from n = —2 to n = —1 for 
the range of different masses mentioned in the figure. In 
other words, the 'bundle' of curves corresponding to each 
value of a{M) covers a large range of M, n and <jg. The 
fact that the curves fall into tight bundles defined only 
by o~{M) confirms what we have inferred from Eq. 1: for 
a given value of f2 , the redshift evolution is almost com- 
pletely independent of n and and does not depend directly 
(78 . On the other hand, there is a significant difference be- 
tween the two cosmological models - as much between the 
open and critical models shown as between a m = 0.66 and 
<r m = 1.0 — a difference significant enough to potentially 
discriminate between the two cosmologies. 
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3. Comparison with observations 

In order to apply this technique, it is important to notice 
that the mass in the PS formula corresponds to a fixed 
contrast density and, therefore, represents very different 
objects at different redshifts. As this mass is not directly 
observable, we must resort to some other, more observable 
cluster quantity. To this end, we use cluster temperature 
and introduce an evolution coefficient: 
n{> T,z) 



c T (z) 



(4) 



n(> T,z = 0.) 

which we will consider for two temperatures - 4 and 6 keV. 
One must then take into account the fact that clusters 
with identical temperatures at different redshifts corre- 
spond to different masses (in the PS language - see Eq. 
1); thus, the evolution expressed in terms of temperature 
could in principle be sensitive to the spectrum. 

We estimate our modeling uncertainty using the re- 
sults of Oukbir et al. (1997, hereafter OBB) and OB. For 
f2o = 1, we allow a(M — 1O 15 M ) to cover the range 
0.55 — 0.65, i.e. a% in the range 0.53 — 0.625 in agreement 
with Viana and Liddle (1996). Rather than the best fit- 
ting value of n = —1.8 given by OBB, we use n — —1.4, 
because it is closer to a F-CDM model with T = 0.25; 
this reduces the amount of evolution by a factor of 2 at 
z ~ 0.5, relative to the n = —1.8 case. For the open model, 
we set f2 = 0.3 and consider two extreme cases with 
a(M = 10 15 M Q ) = 0.78 (according to the results of OB) 
: n = -1.8 with cr 8 = 1.07, and n = -1.2 with ct 8 = 0.94. 
For these ranges of parameters, we examine the redshift 
evolution of the cluster number density for 4 keV and 6 
keV clusters. The corresponding range of predictions for 
ct are presented as the grey areas in Fig. 2. Notice that 
(t(10 15 M q ) differs slightly between the two models (see 
the previous section), increasing the evolutionary differ- 
ence between them. As one can see, the uncertainty for 
£lo = 1 is rather large, but the probe can certainly dis- 
criminate between a low-density and a high-density uni- 
verse. 

It is difficult to directly apply this test to present-day 
X-ray cluster samples, because this requires knowledge 
of the temperature distribution function at high z. The 
only well controlled high-redshift sample of X-ray clusters 
is the EMSS (Gioia & Luppino 1994). It has been stud- 
ied and modeled in detail by OB. They concluded that, 
in order to self-consistently model X-ray clusters in an 
open universe, one must introduce negative evolution in 
the temperature-luminosity relation (i.e., at a given tem- 
perature, clusters are less luminous in the past). The rea- 
son is that the EMSS sample provides definitive evidence 
for negative evolution of the X-ray luminosity function 
(see the following discussion) , while an open cosmological 
model would predict an X-ray temperature function with 
little evolution. 

Recently, several authors have quoted numbers for the 
redshift evolution of the cluster number density. Carlberg 
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Fig. 2. Relative evolution of the abundance of clusters 
above a given temperature. The continuous and dashed 
lines show the cluster abundance evolution for T > 4 keV 
and T > 6 keV, respectively The grey area delimits our es- 
timate of the modeling uncertainty, taken from Oukbir & 
Blanchard (1997) and Oukbir et al. (1997). The triangles 
are from the observations as given by Carlberg al. (1997) 
and Bahcall ct al. (1997). The point at z = 0.66 is derived 
from the luminosity function of the EMSS. The open tri- 
angles correspond to clusters with T > 6 keV, while the 
filled triangle is for T > 4 keV, assuming no evolution in 
the L x — T x relation. 

et al. (1997) have estimated the number density of CNOC 
clusters with velocity dispersions > 800 km/s. They find 
n(z = 0.22) = 4.38 x 10" 8 - and n(z = 0.45) = 1.13 x 10" 8 -. 
We may convert the velocity dispersion to an X-ray tem- 
perature of T x sa 5. keV (in agreement with their luminos- 
ity of £[o.3-3.5] ~ 4.10 erg/s) using the conversion pro- 
vided by Sadat et al. (1998), which shows good agreement 
with recent ASCA measurements (although a few clusters 
appear discrepant). Henry (1997) provides the first actual 
estimate of evolution of the temperature distribution func- 
tion, although at moderate redshift (z w 0.35); the data 
seem to indicate a significant amount of evolution. Fan et 
al. (1997), using the CNOC sample, find 



n(z = 0.5) 

n(z = 0.) 



0.2 



(5) 



for clusters of mass M1.5 = 6.310 1 h _1 M Q within a phys- 
ical radius of 1.5h _1 Mpc. This corresponds to an approx- 
imate temperature of 4.5 keV for a virialized cluster (in- 
dependent of redshift). 

The abundance of X-ray clusters at redshift 0.66 can 
be estimated from the EMSS (Luppino and Gioia 1995): 



N(L [0 . 



3-3.5] 



> 5.10 4 



1.10 



-s 



(6) 



In the absence of evolution in the L x — T x relation , such 
clusters would have temperatures greater than 5 keV (Ar- 
naud and Evrard 1997). The abundance of clusters de- 
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duced from the temperature distribution function at z = 
is, rather surprisingly, highly uncertain (see, for instance, 
Table 1 in Carlberg et al., 1997). To lower this uncertainty, 
we estimate the present-day abundance of similar clusters 
from the BCS luminosity function, which is constructed 
from a much larger cluster sample, (Ebeling et al., 1997). 
In the ROSAT band - [0.1, 2.4] keV - 4 to 6 keV clusters 
have a luminosity greater than 4.10 44 erg/s/cm 2 , yielding 
giving: 



N(> L) ~ 0.6 10~ 7 



ut(z = 0.66) 
n T (z = 0.00) 



0.16 ±0.06 



(7) 



which is direct and clean evidence for some kind of evo- 
lution. The above density will serve as our reference for 
the abundance at z = for the CNOC clusters: ct(~z = 



+0.075 
045' 



0.27) PS O.Utoil and c T (z = 0.45) ps O.IIIq 

As we have already mentioned, open models (f2 ~ 0.2), 
for which the temperature distribution function shows lit- 
tle evolution, cannot be consistent with the EMSS dis- 
tribution unless there is strong negative evolution of the 
luminosity-temperature relation: whatever the value of 
f2o, the properties of the cluster population (cither the 
number density or the luminosity-temperature relation) 
must evolve in order to explain the EMSS redshift dis- 
tribution. One may wonder whether a bias in the EMSS 
sample could lead to a severe underestimation of the clus- 
ter abundance at large z. This seems rather unlikely, for 
at such redshifts clusters are almost point-like compared 
to the size of the detection cell (5'); furthermore, no sys- 
tematic bias has been found in the photometry (Nichol et 
al., 1997). 

These numbers already give interesting insight con- 
cerning the density parameter of the universe. It is clear 
from Fig. 2 that the critical model is favored over a low- 
density model, according to the cluster abundances re- 
ported in the recent literature; however, a note of cau- 
tion: it must be remembered that in all cases, the data 
were analyzed assuming, either implicitly or explicitly, a 
non-evolving relation between temperature and luminos- 
ity. It is for this same reason that our present conclu- 
sions are exactly the same as those given by Oukbir and 
Blanchard (1997): under the assumption of a non-evolving 
temperature-luminosity relation, the EMSS redshift dis- 
tribution of X-ray clusters favors a high density universe. 
This result is supported by the additional information that 
available data on distant X-ray clusters does not demon- 
strate any sign of the strong negative evolution of the 
luminosity-temperature relation needed to save the open 
model (Sadat et al., 1998) (this is independent of the pos- 
sible addition of a cosmological constant). This additional 
piece of information is critical to the conclusion, because 
without it, we have no way of understanding the flux lim- 
ited selection of the EMSS in terms of temperature. 



4. Conclusion and discussion 

The purpose of this letter was to clarify the nature of the 
evolution of the cluster temperature distribution function. 
As we have seen, this evolution depends primarly on the 
amplitude of the fluctuations on the scale under consider- 
ation, <r(M), and the cosmological background, f^o- This 
is the origin of the robustness of the cosmological test 
initially proposed by Oukbir and Blanchard (1992). The 
EMSS redshift distribution, as modeled by Oukbir and 
Blanchard (1997), combined with the absence of observed 
negative evolution in the temperature-luminosity relation 
provided the first evidence for a high density universe from 
this technique (Sadat et al 1998). Our analysis leads to a 
similarly high value for the density of the universe. During 
the submission of this letter, we learned that similar con- 
clusions were reached by two other groups who included 
ROSAT cluster redshift distributions (Borgani et al., 1998; 
Reichart et al, 1998). Because this test is primary sensi- 
tive to the dynamical behavior of the universe as a whole 
(through the growth rate of linear density fluctuations), 
we consider this to be the strongest evidence in favor of a 
critical density universe presently available. 
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